Persistent currents in one dimension: the other side of Leggett's theorem 
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We discuss the sign of the persistent current of A'' electrons in one dimensional rings. Using 
a topology argument, we establish lower bounds for the free energy in the presence of arbitrary 
electron-electron interactions and external potentials. Those bounds are the counterparts of up- 
per bounds derived by Leggett. Rings with odd (even) numbers of polarized electrons are always 
diamagnetic (paramagnetic). We show that unpolarized electrons with A'^ being a multiple of four 
exhibit either paramagnetic behavior or a superconductor-like current-phase relation. 
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A persistent current fH flows at low temperatures in 
small conducting rings when they are threaded by a mag- 
netic flux (j). This current is a thermodynamic effect 
which is deeply connected to the presence of quantum co- 
herence. Its magnitude can be expressed as / = —dF/d(j) 
in terms of the free energy F. Persistent currents have 
been in the focus of an intensive theoretical activity (see 
for example SBflSSi)- Nevertheless, the understand- 
ing of the exp erimentally measured currents in metal- 
lic 0, H, 0, and semiconductor 11 , l3| rings remains 
incomplete. In particular, the observed currents in dif- 
fusive rings are an order of magnitude larger than the 
value obtained from theories neglecting electron-electron 
interactions. 

The persistent current itself generates a magnetic fleld 
which is detected in the experiments. Besides the mag- 
nitude of the current, the sign of this magnetic response 
is of particular interest. In the absence of interactions, 
mesoscopic sample-to-sample fluctuations between para- 
magnetic {F{(t)) < F{0)) and diamagnetic (i^((/)) > F{0)) 
behavior are expected. Taking into account repulsive (at- 
tractive) interactions leads to the prediction of a param- 
agnetic (diamagnetic) average response in ensembles of 
diffusive rings. However, only diamagnetic signals have 
so far been observed in experiments. 

A very general theoretical result in this domain is the 
theorem by Leggett [l3| which states that the sign of the 
zero-temperature persistent current for spinless fermions 
(fully polarized electrons) in one dimensional (ID) rings 
is given by the parity of the particle number. This result 
appears in the form of upper bounds for the ground state 
energy. Leggett's theorem holds for arbitrary potential 
landscape and electron-electron interactions. 

In this letter, we perform two tasks. Firstly, we estab- 
lish general lower bounds of the free energy that can be 
seen as the counterparts of Leggett's upper bounds. Our 
result is valid at any temperature and can be general- 
ized to include the spin degree of freedom. Secondly, we 



use the density matrix renormalization group (DMRG) 
method [3| to study scenarios allowed by those bounds. 
We start with the Hamiltonian 
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describing A^ spinless electrons in a ID ring of size L, 
where pi is the momentum of electron i, A(r) the vec- 
tor potential corresponding to a magnetic flux 0, V{r) 
an external potential and U{r) the electron-electron in- 
teraction. A gauge transformation allows to replace the 
vector potential by flux-dependent boundary conditions 

^'(L,r2,...,rjv) = e'**(0,r2,...,rAr) 

with <I> = 2Tr(j)e/h for the antisymmetric wavefunction \E'. 
Leggett's argument 13| is based on the ansatz 



^v{ri, ■ . . ,r7v) = e 



for the variational ground state at (f> ^ in terms of 
the exact ground state wavefunction at $ = 0. The 
energy corresponding to this ansatz is 



= i;($ = 0) 
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where E{^) is the exact ground state energy. The phase 
is a symmetric function of its arguments and should go 
from to $ as ri goes around the ring from to L. In 
this loop, ri crosses the coordinates of the A^ — 1 other 
particles. At each crossing the sign of ^'o changes. Since 
^'o must be back to its initial value at the end of the loop, 
for even N there must be at least one more sign change. 
The function can be chosen such that its gradient is 
concentrated around this additional point where van- 
ishes, yielding Ey — E{0). This variational energy pro- 
vides an upper bound E{^) < E{0) for the ground state 
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FIG. 1: Four different paths R{r) of the path integral ([T]). 
The winding numbers m are equal to (a), (b), —2 (c) and 
+1 (d). Path (b) does not contribute in ID rings because two 
particle lines cross, which is forbidden by antisymmetry. 



energy and the system is paramagnetic. A similar argu- 
ment leads to < E{7r) for odd N. In the sequel, we 
show that these upper bounds can be complemented by 
corresponding lower bounds. 

We consider the partition function Z = TrAe~^^ 
where Tia is the trace over antisymmetrized many-body 
states. One can divide the inverse temperature (3 'm K 
increments Ar — j3/K and insert 1 = J dR\R){R\ at 
each step, where \R) = |ri, r2, . . . , r^) denotes the list of 
positions of all particles. The result reads 

P Ri,...,Rk 

■ • • [e \Rk-iRk [e \r.kP(Ri) : 

where P is a permutation of the N particles and \P\ = ±1 
its signature. This path integral formula can be formally 
rewritten in continuous form as 



Z = J2\P\ [ VR{T)e- 

ID J 



S\R{r-)\ 
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where S\R{t)\ is the action of the path R{t) which satis- 
fies i?(/3) = P[i?(0)]. Examples for such paths are shown 
in Fig. [T] The effect of a magnetic flux on the path 
integral is a phase ±$ each time a particle crosses the 
boundary, yielding 
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where the winding number TO[i?(T)] counts the total num- 
ber of boundary crossings contained in R{t) (+1 from left 
to right and —1 in the opposite direction). With F ~ 
-\ogZ/p, Eq. (ED leads to d^F/d^\^^ = {m^[R])/[3 



which is widely used to determine the superfluid fraction 
in bosonic systems, in particular in the context of quan- 
tum Monte-Carlo simulations Anyway, Eq. ([2]) can 
be expressed as 
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with the partial partition functions Z,„ for a given m. 

In ID rings, most of the permutations P do not con- 
tribute to the path integral since they would correspond 
to paths where two particles cross at some point. Those 
paths (see Fig. [T})) have zero contribution due to anti- 
symmetry. As a result, only paths with cyclic permu- 
tations contribute to the path integral. There is only 
one single cyclic permutation Pm per winding number m 
(see Fig. [1]), whose signature determines the sign of the 
corresponding Z„i. Since the signature of a cyclic per- 
mutation of N particles with winding number m is given 

by = (_l)m(Ar-m)^ Yifwe 
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If is odd, then ra(N — m) is even and is positive 
for all m. Since Z^n = Z-^^ we have 

OQ OC 

Z($) = Zo+2 E ^mCOs(m$) < Zo+2 E \Z,n\ = Z{0) 

m— 1 m—1 

SO that F{^) > F{0), and the system is diamagnctic. 
When N is even, m{N — m) has the same parity as m 
and Zjn = (— leading to 

OC 

Z($) < Z0 + 2E \Z.n^\=Z{^) 
m—1 

such that F{^) > F{tt). At zero temperature, the above 
inequalities together with those obtained by Leggett read 

E{0) < E{^) < E{-k) for N odd, (3) 
E{tt) < E{^) < E{0) for N even. (4) 

The previous results can be partially extended to un- 
polarized ID systems with N-^ up and down electrons 
under a spin conserving Hamiltonian. Introducing wind- 
ing numbers m| and for the two species, we have 

Z(<i>)= E ^^T.^ie^*'"''""'^- 

Using similar arguments as above, we find Zm^,mi = 
(-l^l-P^Tl+l-^^illZm^^mJ and conclude that 

F{0) < F($) for and Ni odd, (5) 
F{tt) < P($) for and Ni even. (6) 

Those lower bounds are valid at arbitrary temperature. 
Eqs. ^ and ^ for the polarized case and ^ and ^ for 
unpolarized electrons are the central results of this work. 
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FIG. 2: Two possible flux-dependences of the free energy F 
for N = 4n with even 7V| = A*'|, corresponding to paramag- 
netic (i) and diamagnetic (ii) behavior. 

We now consider rings with a number of electrons 
N = Ni + Ni = An being a multiple of four, and char- 
acterize the flux-dependence of the free energy using the 
curvatures 

Co(7r) = 9^-^/9$^ I $=o(Tr) 

at $ = and tt. In the unpolarized case {N^ = — 2n) 
Eq. ^ applies, the free energy assumes a minimum at 
$ = TT, and Ctt is positive. The curvature at zero flux Cq 
is not constrained such that the two scenarios illustrated 
in Fig. [2] remain, (i) If Co is negative, the system is 
paramagnetic and lowers its energy in the presence of a 
magnetic flux, (ii) If Co > 0, the system is diamagnetic. 
In this case the usual sinusoidal current-phase relation 
F oc — cos$ is prohibited by > 0. There is a strong 
contribution of the second harmonic in $ and it exists at 
least one flux value < $c < tt where the free energy is 
maximum and the persistent current vanishes. Without 
interaction, one finds case (i). Case (ii) is characteristic, 
for instance, of superconducting fluctuations induced by 
a weak attractive interaction (in a superconductor, F oc 
-cos(2$) and $c = 7r/2). 

The overall conclusion is that a simple diamagnetic 
response is prohibited. The system must either be 
in a paramagnetic state or have a superconducting-like 
current-phase relation. 

In order to illustrate the behavior of the curvatures Cq 
and Ctt, we consider ID tight-binding rings of M sites, 
described by the Haniiltonian 

M 

a j=l 
M 

+ (^o^^J^T"J,i + UiUj+inj) . (7) 

Here, ^ creates an electron with spin a = {|, J,} on site 
j, the random on-site energies Vj,^ are drawn indepen- 
dently from the interval [— 1^/2, and the kinetic 
energy scale is given by the hopping amplitude t. We 
include Hubbard on-site and nearest-neighbour interac- 
tions of strength Uq and Ui, respectively. The density op- 
erators are defined as nj^cr = Cj ^Cj ^ and rij = rij^-^ 
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FIG. 3: Curvatures Co at zero temperature for disordered 
quarter-fiUed Hubbard rings of different sizes M — 2N , as 
a function of Uo, a.t U\ = Q and A^^ = -^i, for one disorder 
realization with W — t. 

and the boundary condition cm+i,(t — e**ci_cr accounts 
for a magnetic flux threading the ring. 

We use the DMRG algorithm adapted to disor- 
dered systems [31 to calculate the ground state energies 
and the zero-temperature curvatures 17] for $ = and 
$ = TT, fully taking into account the many-body corre- 
lations. For the largest systems 750 states per block are 
kept in the DMRG iterations. All of the numerical results 
should, and do, satisfy the relations © and (O. 

In Fig.[3l we show the effect of an attractive interaction 
on the curvature Co, for single realizations of disordered 
quarter-filled Hubbard rings of sizes up to M = 2N = 32 
sites with W — t and Ui = 0. For even = Ni 
{N ~ An), an attractive interaction reverses the sign 
of Co. Hence the interactions induce a transition from 
scenario (i) to (ii) of Fig. [2] and drive the system from 
paramagnetic towards "superconducting". In contrast, 
the sign of Co for odd N-^ = Ni {N = An + 2) remains 
positive for all values of Uq, as dictated by ([5]). 

A molecular realization of the N ^ An case is cyclooc- 
tatetraene (CsHg, see the inset of Fig. [3]), which consists 
of a ring of eight carbon atoms with eight tt electrons. In 
Fig. m we plot the curvatures Cq and Cjr as a function 
of the nearest-neighbour interaction strength Ui for the 
model Hamiltonian ([7]) with M — 8 and W = 0, using the 
parameters t = 2.64 eV and Uq = 8.9 eV given in Ref. [31 
for cyclooctatetraene. Depending on the strength of Ui 
(and the neglected longer range parts of the interaction) , 
the system can undergo a transition [Toj from a para- 
magnetic spin-density- wave {Ui < 4.6 eV and Co < 0) 
to a diamagnetic charge-density-wave {Ui > 4.6 eV and 
Cq > 0) groundstate. 

The pres ence of paramagnetism or orbital ferromag- 
netisrn [ 2 p| in cyclooctatetraene is still a matter of de- 
bate [iSl, l2l|. A ferromagnetic instability can occur in 
small paramagnetic rings provided their inductance is 
large enough. Then, a magnetic field fluctuation gen- 
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FIG. 4: Curvatures Co and Ctt as a function of Ui for an 
unpolarized clean ring with M = 8 at zero temperature, for 
the parameters of cyclooctatetraene (A^ — 8, t = 2.64 eV, 
Uq = 8.9 eV). Inset: Sketch of the C atom configuration in 
cyclooctatetraene. 

erates a current which reinforces the magnetic field. For 
small flux, the persistent current is / « Coe$/fi. and the 
magnetic flux generated by this current is « noeLI/h 
(L » 1 nm is the typical circumference of the molecule 
and hqL is its typical inductance). The instability occurs 
when the magnetic field generated current amplifies field 
fluctuations, i.e. when Cq < and X — fiQ\CQ\e'^ L / h? > 
1. In our model, the most pronounced negative curva- 
ture Co ~ — 5eV occurs in the spin-density- wave regime 
close to the transition at J7i w 4.6 eV. We therefore have 
X < 3 X 10~^, far from the ferromagnetic instability. 
Moreover, Jahn- Teller distortions in cyclotetraene have 
been predicted j2li | which would reduce the curvatures. 

An alternative to small molecules are rings made in 
semiconductor heterostructures. Important progress has 
been realized in the fabrication of such systems and ID 
rings with a single conduction channel could be produced. 
Estimates show that one still has X 1. Therefore, 
no orbital ferromagnetism is to be expected, reminiscent 
with the result for large ID Luttinger liquid rings [22]. 

However, as illustrated by the results shown in Fig. [3l 
we predict that changing the number of electrons by two 
(with a back gate for instance) leads to a change of either 
the sign of the magnetic response (paramagnetism) or 
the periodicity of the response to a magnetic field ( "su- 
perconducting like"). This could allow to detect small 
attractive interactions in those systems, whose presence 
is suggested by the fact that so far only diamagnetic re- 
sponses have been observed in multichannel rings. 

In conclusion, we have established general relations for 
the sign of persistent currents in ID systems that can 
be considered as the counterpart of Leggett's theorem 
for spinless fermions. Our theorem includes spin and is 
valid at arbitrary temperature. For spinless fermions, 
our relations imply that interactions and disorder cannot 
affect the sign of the persistent current. In particular. 



the possibility discussed by Leggett of having maxima of 
the energy at both, integer and half integer flux quan- 
tum, corresponding to a paramagnetic signal in an as- 
sembly of rings, is ruled out by ([3|) and (jl]). For electrons 
with spin, when Leggett's theorem does not apply, only 
the lower bounds ([5]) and ^ for the free energy remain. 
We showed that this allows for a "superconducting-like" 
current-phase relation when attractive electron-electron 
interactions change the sign of the unconstrained cur- 
vature. Our topological argument provides a rigorous 
justification for the phenomenological Hiickel rule which 
states that cyclic molecules with 4n + 2 electrons like 
benzene are aromatic while those with An electrons are 
not. 
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